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ABSTRACT 

The recent measurements of weak lensing by large-scale structure present significant 
new opportunities for studies of the matter distribution in the universe. Here, we 
present a new cosmic shear survey carried out with the Echelle Spectrograph and 
Imager on the Keck II telescope. This covers a total of 0.6 square degrees in 173 
fields probing independent lines of sight, hence minimising the impact of sample vari- 
ance. We also extend our measurements of cosmic shear with the William Herschel 
Telescope (Bacon, Refregier & Ellis 2000) to a survey area of 1 square degree. The 
joint measurements with two independent telescopes allow us to assess the impact 
of instrument-specific systematics, one of the major difficulties in cosmic shear mea- 
surements. For both surveys, we carefully account for effects such as smearing by the 
point spread function and shearing due to telescope optics. We find negligible resid- 
uals in both cases and recover mutually consistent cosmic shear signals, significant 
at the 5.1(7 level. We present a simple method to compute the statistical error in 
the shear correlation function, including non-gaussian sample variance and the co- 
variance between different angular bins. We measure shear correlation functions for 
all fields and use these to ascertain the amplitude of the matter power spectrum, 
finding (78 (^)°'^^ = 0.97 ± 0.13 with 0.14 < < 0.65 in a ACDM model with 
r = 0.21. These 68% CL uncertainties include sample variance, statistical noise, red- 
shift uncertainty, and the error in the shear measurement method. The results from 
our two independent surveys are both consistent with measurements of cosmic shear 
from other groups. We discuss how our results compare with current normalisation 
from cluster abundance. 

Key words: cosmology: observations - gravitational lensing, large-scale structure of 
Universe. 



1 INTRODUCTION 

Weak gravitational lensing of background galaxies by inter- 
vening large-scale structure ('cosmic shear') provides direct 
information about the total mass distribution in the uni- 
verse, regardless of its nature and state. A measurement of 
cosmic shear thus bridges the gap between theory, which 
is primarily concerned with dark matter, and observation, 
which generally probes only luminous matter. The recent 
measurements of coherent distortion of faint galaxies by sev- 
eral groups (van Waerbeke et al. 2000; Bacon, Refregier & 
Ellis 2000 [BRE]; Wittman et al. 2000; Kaiser, Wilson & 
Luppino 2000; MaoU et al. 2001; Rhodes, Refregier & Groth 
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2001; van Waerbeke et al 2001; Hammerle et al. 2001; Hoek- 
stra et al. 2002; Refregier, Rhodes & Groth 2002) has trig- 
gered great interest in the provision of new constraints on 
the amount and distribution of dark matter, together with 
measurements of several cosmological parameters. 

If intrinsic galaxy orientations are essentially random 
in a given survey (which requires the survey to be suffi- 
ciently deep; see Brown et al. 2000; Catelan et al. 2000; 
Heavens et al 2001; Croft & Metzler 2001; Crittenden et 
al 2001), any coherent alignment must arise from distor- 
tion due to weak lensing. Light paths from galaxies pro- 
jected close together on the sky pass through, and are grav- 
itationally distorted by, the same dark matter concentra- 
tions. This coherent distortion contains valuable cosmolog- 
ical information (eg. Bernardeau et al. 1997; Jain & Seljak 
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1997; Kamionkowski et al. 1997; Kaiser 1998; Hu & Tegmark 
1998). In particular, the variance of the distortion field mea- 
sures the amplitude of density fluctuations (~ as,^^). This 
shear measurement is free from assumptions about gaussian- 
ity or the M — T relation, and whilst the shear-based mea- 
surement is currently comparable in precision to that from 
local cluster abundance, further progress is limited solely by 
the number of fields observed. 

The validity of results from cosmic shear surveys de- 
pends sensitively on the treatment of systematic errors. A 
further issue arises from sample (or 'cosmic') variance, the 
impact of which can be limited by using numerous indepen- 
dent sightlines to complement panoramic imaging of a few 
selected areas. With these motivations in mind, we present 
a comparison of the cosmic shear observed with two inde- 
pendent instruments (Keck and WHT), using two different 
survey strategies. 

In this paper, we describe the first cosmic shear mea- 
surements with the 8' x 2' Echelle Spectrograph and Imager 
(ESI) on the Keck II telescope. This Keck survey reaches a 
depth of z ~ 1.0, comparable to other recent cosmic shear 
surveys (e.g. van Waerbeke et al 2001, Bacon et al 2000). 
However, the much faster acquisition of fields with ESI in 
comparison to 4m telescopes such as the William Herschel 
Telescope allows us to obtain very many more fields (173 in 
the final survey) . This improves the cosmic shear signal mea- 
sured, by minimising the contribution to noise of the sample 
variance, i.e. the error upon the mean lensing signal result- 
ing from the measurement of shear on only a limited number 
of fines of sight. In order to measure the cosmic shear signal, 
we analyse the correlation function of the shear on various 
scales, and obtain constraints on cosmological parameters 
using a fit to theoretical predictions upon varying these 
parameters. 

In addition to this investigation, we extend our original 
detection of cosmic shear on the 4.2m WHT (in BRE), to 
a measurement of the correlation function of the distortion 
field for this dataset. A more precise measurement is afforded 
by the increase in the number of WHT fields to a total of 
20, with a further increase in area due to the larger 16' x 16' 
size of field with the new WHT mosaic camera. 

This paper is organised as follows. In §^ we discuss our 
observational strategy for measuring cosmic shear with both 
telescopes. We then describe (§^ the procedure for observ- 
ing with the Keck telescope, and explain the method used 
for reducing the data. We follow this by a treatment of the 
systematic effects associated with the Keck data, such as 
shear distortion from the camera and the anisotropic smear 
due to telescope tracking. 

Having dealt with the Keck data up to the stage of mea- 
suring galaxy shapes, we review the procedure for obtaining 
the WHT data in M followed by the approach adopted for 
data reduction and removal of systematic effects. 

We proceed to develop the correlation function formal- 
ism in §pl which we use for precision measurements of the 
cosmic shear, in §^ We then interpret the shear signal in 
terms of cosmological models in order to obtain limits on 
the cosmological parameters as and Qm- Our conclusions 
are summarised in sM 



2 SURVEY STRATEGY 

The aim of our Keck and WHT surveys has been to acquire 
deep {z ~ 1) fields representing numerous independent lines 
of sight, sufficiently scattered to sample independent struc- 
tures and thus minimise uncertainties due to sample vari- 
ance. These lines of sight must be chosen in a quasi-random 
fashion, without regard to the presence or absence of mass 
concentrations, in order to obtain a representative sample 
of the mass fluctuations in the universe. Here we describe 
the strategy adopted for the two surveys, based on that in 
BRE. 

The survey fields for both Keck and WHT were selected 
by choosing a sparse (> 2° separation for statistical inde- 
pendence) grid of coordinates spanning the range accessible 
to the telescope at a given time. We tuned the Galactic lati- 
tude of the grid to ensure ~50 unsaturated stars within the 
Keck field of view, and ~200 in the WHT fields, so that the 
anisotropic PSF and the camera distortion could be care- 
fully mapped (see below). The STScI Digitised Sky Survey 
was then used to find an appropriate final field near each 
set of coordinates, avoiding stars brighter than i? < 11 in 
the APM and GCC catalogues, to prevent large areas of 
saturation or ghost images. 

As a final constraint, each field was observed within 20° 
of zenith for both telescopes. This minimises smearing due to 
atmospheric refraction for Keck, which does not have an At- 
mospheric Dispersion Corrector (see section ^ and figure |l^ 
for confirmation that this is not a limiting systematic). For 
both telescopes the constraint minimises any image distor- 
tion associated with telescope or instrument fiexure. Figure 
|l] shows the positions of the resulting selected survey fields 
on the sky. 

We must now determine the depth to which to observe 
these fields. We have shown in BRE that the cosmic shear 
signal is measurable with WHT images having a 1 hour ex- 
posure length, corresponding to a usable galaxy sample with 
median source redshift Zs — 0.8. We have further demon- 
strated in Bacon et al (2001) that longer exposures do not 
improve the signal greatly, as beyond this depth galaxy 
shapes are seriously degraded by a typical ground-based 
PSF. On the other hand, shorter exposures face the dan- 
ger of considerable contamination by intrinsic alignments 
of galaxies (e.g. Heavens et al 2001). We therefore aimed 
to probe the same redshift range with our new Keck sur- 
vey. An exposure time of 10 minutes was calculated to allow 
5(7 detections of point sources at i? = 26, given the optics 
throughput, sky background in the Keck R filter frequency 
range and the quantum efficiency of the ESI CCD. How- 
ever, as we shall see, the better seeing during our observa- 
tions with Keck results in a slightly fainter magnitude limit 
than the WHT survey; this is entirely acceptable, as we can 
compare results by scaling the predictions according to the 
equation oc z"'^ (see e.g. BRE). 

The fact that gravitational lensing is achromatic per- 
mits us a free choice of photometric band for our observa- 
tions. However, we should note that R and / afford the most 
efficient deep imaging in a given exposure time. Due to fring- 
ing in the / band with the EEV CCD in use at the WHT, 
we choose to image in R with both telescopes. The Keck R- 
band is a specially constructed filter with similar throughput 
and spectral range to the WHT Harris R; the slightly dif- 
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- ■ Keck fields 
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Figure 1. The sky location of tiie cosmic shear fields. Galactic 
latitudes of 0°, ±20°, ±40°, ±60° are shown as contours; the 
Galactic centre and poles are shown as a cross. 




Figure 2. Histogram of seeing (FWHM) in all survey shear fields 
after stacking dithered exposures. Keck data are shown in white; 
WHT data are shown as shaded. Each 8' X 16' WHT chip is 
counted separately. Combining survey area from the two tele- 
scopes, the median seeing is 0.73"; no data with seeing worse 
than 1" is used. 



Keek fields 
(inediyii fl=24.4) 



WHT fields 
(median fl=23,e 



Median dcpfh of oalaxy cafalogue {R mags) 

Figure 3. Median reddening-corrected i?-band magnitude of the 
galaxies in each field that are used in the final object catalogue. 
The median depth of the combined surveys is _R=24.0. 



ferent galaxy distribution probed by this filter is taken into 
account later in our redshift uncertainty estimates. 



3 KECK DATA 
3.1 Observations 

As we briefly mentioned above, the great advantage which 
Keclc presents in measuring cosmic sliear is its speed in 
achieving the necessary depth; 10 minute exposures reach 
i? = 26 (contrast the 1 hour necessary for such a depth 
with WHT). We can therefore obtain very many indepen- 
dent lines of sight per night, thus reducing the impact of 
sample variance. 

We observed 173 2' x 8' ESI fields using a specially made 
R filter (A = 6657A, eflfective FWHM(A) = 1200A), over 
the course of 6 nights in June 2000, November 2000 and 
May 2001. The necessary imaging observations were done 
as a good seeing override on an independent spectroscopic 
programme. The pixel size of 0.153" for ESI is considerably 
finer than that for WHT, allowing better sampling of the 
galaxy images. 

For each of our fields, three separate exposures were 
taken, each offset by 5". This enables the continual re- 
calibration of optic al distortions in the telescope and cam- 
era (see section 3.3), and the removal of cosmetic defects and 



cosmic rays. All the fields were observed as they passed near 
the meridian (in order to minimise image distortion from the 
atmosphere and from telescope or instrument flexure) but 
no closer than 5° to zenith (to minimise sky rotation and po- 
tential tracking errors on this alt-az telescope). Bias frames, 
sky flats and dome flats were acquired at the start and end 
of each night, and standard star observations were taken 
regularly throughout each night. The telescope was focused 
several times per night to minimise camera distortions. 

As exposures were observed, we used our real-time 
SExtractor-based (Bertin & Arnouts 1996) software on re- 
cently completed exposures in order to monitor PSF size 
and stellar ellipticity. Figure ^ records the seeing values for 
all fields; we found that the median seeing for the observa- 
tions was 0.59", with 75% in seeing better than 0.7". As we 
discuss below, this excellent quality data yields a low level 
of noise on the estimates of cosmic shear. 

Furthermore, the rms stellar ellipticity for the images 
was Get ~ 0.035. This relatively low value facilitates the 
necessary PSF corrections (see e.g. Bacon et al 2001, Er- 
ben et al 2001). The stellar ellipticity on Keck fields is often 
found to be due to tracking errors, seen as a uniform PSF 
anisotropy across the field of view (see figure |^ for an ex- 
ample and discussion in §3.4). In exposures with excellent 



tracking, the dominant effect appears to be astigmatism due 
to the fact that the CCD is slightly skewed with respect 
to the focal plane, thus probing optical conditions slightly 
above and below the focus. 

In terms of depth, with magnitudes calibrated using 
photometric standards, we found that the median magni- 
tude of all galaxies detected was i?=25.1, with an imcat 
signal-to-noise of 5.0 being reached for galaxies at J?=25.8. 
We keep ~ 27.5 resolved galaxies per squ are arcminute in 
our final object catalogue (see section 



3.4 for the selection 



criteria), which corresponds to a final magnitude median of 
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R ~ 24.4 ± 0.2, including reddening corrections. The distri- 
bution of magnitudes for the final galaxy catalogue extends 
substantially fainter than this median, with the galaxy count 
as a function of magnitude dropping to 50% of maximum at 
R = 25.2. According to e.g. Cohen et al (2000), our median 
magnitude of 7? = 24.4 corresponds to a median source red- 
shift of Zs ~ 1.0. We can estimate the error on this by tak- 
ing Cohen et al's measured interquartile range of redshifts 
for galaxies at this magnitude, /Q=l. 6-0. 55=1. 05, and us- 
ing the estimate of the error on the mean for a gaussian 
distribution, a — IQ/{1.35^/n) = 0.1, where the number of 
galaxies with redshifts measured n — 49. Thus we will use an 
estimate of the median redshift Zs — 1.0 ±0.1. This quantity 
may be subject to additional sample variance, which would 
increase the error; we will therefore quote our redshift er- 
ror separately to other statistical errors for comparison with 
future redshift surveys. 



-O.IZ Shear 
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3.2 Data Reduction 

These deep images were reduced following standard meth- 
ods. Bias subtraction of the science exposures and sky-flats 
was individually calibrated for each image using the overscan 
regions at the edge of the chip. The science exposures were 
then trimmed and divided by a median composite (with So- 
clipping) of the flat flelds for a given night. To eliminate any 
remaining background gradients, the science images were 
also divided by a stack of all the (median- normalised) ex- 
posures from that night. In contrast to the WHT data, no 
fringing is observed on the ESI images, due to a thicker CCD 
than that for the WHT camera; this simplifles the reduction 
considerably. 

The multiple exposures for each fleld were aligned by 
cross-matching common objects in SExtractor catalogues 
(Berlin & Arnouts 1996). Typically ~250 objects were found 
in common per Keck fleld, thus providing a stringent es- 
timate of the mean dither offset (to 0.02 pixel accuracy). 
The exposures were then shifted by the required non-integer 
number of pixels via a linear interpolatio n us ing the IRAF 
imshift routine. As we discuss below in 53.3, rotations be- 



tween dithers, and astrometric distortions due to the tele- 
scope and camera optics, are found to be negligible in our 
analysis. 

The resulting registered exposures were then divided 
by their median values to normalise their background levels 
and averaged with Scr clipping to remove cosmic rays and 
cosmetic defects. Remaining problematic areas (edges, bad 
columns, regions containing light leaks from stars outside 
the field of view, spikes from slightly saturated stars, and a 
< 1' square unresponsive region in the corner of the CCD) 
were flagg ed a nd are not used in our cosmic shear analysis 
(see also §3.5). Note that saturation spikes are also largely 
excluded from our object catalogues by an initial e < 0.5 
cut, but a local sky background gradient required that a 
few galaxies within 2" of bright stars also be discarded. For- 
tunately with the small field of ESI we can selectively avoid 
most stars which are bright enough to saturate. After these 
maskings, a reliable ~12.8 square arcminutes per survey field 
remains. 

Having obtained carefully reduced data, the next step 
is to catalogue the shapes of galaxies and to estimate their 
gravitational shear. In order to do this, we must correct for 



0,0 0,5 1,0 1,5 2,0 2,5 3,0 
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Figure 4. Example instrumental shear pattern for ESI. Each 
bar represents the magnitude and orientation of astrometric dis- 
tortions averaged over 20 sets of three dithered exposures. The 
illuminated area of the CCD is rotated from the x — y axes by 7° , 
accounting for the occupied bins seen here. 



any shear introduced by the telescope itself, together with 
tracking errors and atmospheric smearing which convolve 
the galaxy shapes: these effects can mimic a coherent dis- 
tortion which must be removed. We first turn to the issue 
of instrumental optical distortions, mimicking shear. 



3.3 Instrumental Distortions 

Our method for determining the shear fleld induced by the 
telescope and camera, using the offsets in SExtractor cat- 
alogues of our dithered astrometric frames, is fully docu- 
mented in BRE. Figure ^ shows the ESI instrumental shear 
pattern obtained using this method, averaging over 20 flelds 
in order to overcome noise, as the error on the shear is 0.09 in 
a 1' square bin on a given field. We find after this averaging 
that the shear has a mean of 0.2% and is < 0.3% everywhere; 
the shear measured fluctuates by < 0.1% as we average over 
different sets of flelds. This implies that, since our results 
deal with shears of ~ 1%, and since these intrinsic values 
add in quadrature to negligibly affect these shears, we can 
neglect this effect in our analysis. (Note that the ESI fleld 
has an illuminated area which is rotated from the x — y axes 
by 7°, accounting for the slightly slanted geometry on this 
figure.) The magnification and rotation components were 
also found to be < 0.3%, and are consequently negligible 
over such a small field. 
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Figure 5. Example stellar ellipticity pattern for ESI. (Left) Be- 
fore correction; the pattern varies qualitatively between fields 
but mean ellipticity here is typical of the field-to-field average 
of 0.035. (Right) After correction; mean ellipticity reduced to 
<0.001 across the survey, with residual orientations essentially 
randomised. 



3.4 Point Spread Function 

We now turn to more important sources of systematic er- 
ror which need correction: the isotropic smearing of images 
due to the atmosphere and telescope optics, together with 
anisotropic smearing due to tracking errors, dither ahgn- 
ment and optics. 

We used Kaiser's imcat software (Kaiser et al 1995) to 
detect objects on our images, and measure ellipticity, ra- 
dius, magnitude and polarisability (responsiveness to image 
shear/smear; a large galaxy will be less affected by a given 
smear than a small galaxy) for each object. Our procedure 
is described in detail in BRE. 

Noisy detections in our imcat hfindpeaks catalogues 
were removed using criteria as in BRE (r^ > 1.0, v > 15, e < 
0.5). Note that after this, all objects are weighted equally. 
The stars to be used to monitor the PSF were selected from 
the non-saturated locus on a magnitude-radius plane for 
each field. In most observations, the distribution of stellar 
ellipticities over the field is found to be smooth and slowly 
varying (see figure ^for an example). Some exposures taken 
during July 2000 were found to have an unexplained, sharp 
discontinuity in this pattern | of the way up the field. For 
these fields, only objects on one side of the PSF discontinu- 
ity were used. The rms field-to-field stellar ellipticity over all 
used Keck data was found to be Ue* — 0.035; and reduced 
to a negligible <0.001 by our analysis. 

We then measured shear estimators for each galaxy and 
corrected them for convolution with the local PSF as de- 



V 



V 




1.0 1.5 
X (arcmin) 



V 



V 





■1% 

1% 


1% 



















4 



012345678 
y (arcmin) 

Figure 6. Average shear of all galaxies in all Keck fields as a 
function of position on chip. Overall, (71) = —0.02% ± 0.16% 
and (72) = -0.29% ± 0.16%. 



scribed in BRE. To interpolate the PSF to the position of 
each galaxy, we iteratively spatially fitted the measured stel- 
lar ellipticities with a 2-D polynomial, removing extreme 
outliers due to noise and blended images. Similarly, each 
component of the smear polarisability tensor P*^ was fitted 
with a 2-D polynomial. Although individually smooth, the 
PSF patterns were found to have large variations from field 
to field. The degree of fitting polynomial was thus adapted 
to suit each pattern. Typically, a quadratic or cubic compo- 
nent was necessary in the y direction; the much more narrow 
X direction generally required only a constant or linear func- 
tion. We avoided high-order polynomial PSF models which 
diverged towards the edges of the field and would have spuri- 
ously elongated galaxies by over-correcting for the true PSF 
smearing effect (see Massey et al. 2001). 

3.5 Masking of the Field 

As a check for systematic effects associated with data re- 
duction or chip behaviour, we took shear estimators from 
all galaxies in all fields and averaged as a function of CCD 
position (see figure We find that the mean shear for 
our entire ensemble is 71 = —0.0002 ± 0.0016 and 72 = 
—0.0029 ± 0.0016. This is consistent with zero offset in 
the whole ensemble, as we would expect. Figure ^ further 
demonstrates that there is no significant structure in our 
shear values with position on the chip. 

This plot proved very helpful to expose problems with 
the different parts of the analysis (see discussion in Massey 
et al. 2001) . In particula r, th e necessity for stringent mask- 
ing of the field (see also §3.2) quickly became manifest. The 



plot not only mirrored the locations of obvious CCD de- 
fects, but also indicated which other regions of the chip 
were unsuitable for high-precision shape measurements. In 
particular, galaxies at the edges of the data/mask are cut 
apart by the image and appeared aligned with the bound- 
ary. This caused a 2% mean shear offset inside ESI's long 
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and thin field geometry. Other galaxies also appeared spu- 
riously aligned in regions where flat-fielding of any of the 
three dithered exposures was imperfect (e.g. near the edge 
of any one of the dithers; near internal reflections of bright 
stars outside the field of view; or through cirrus). In these 
instances, the dithers have slightly different background lev- 
els and image co-addition leads to residual ellipticities in 
galaxies around these locations. If they are not excluded, 
the mean shear throughout the field is again offset by 
we found that the masking of approximately one quarter of 
the CCD reduced these effects to negligible levels in ~12.8 
square arcminutes of selected data per field. 
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4 WHT DATA 



4.1 Observations 



The WHT data used includes the 13 8' x 16' fields used in 
our previous detection of cosmic shear (BRE), as well as 7 
fields observed in June 2000 with the new larger field of view 
for the Prime Focus Camera (16' x 16'). Thus the total area 
of the combined WHT survey is 1.0 square degree. 

The fields were observed using the Harris R filter with 
the William Herschel Telescope Prime Focus Camera. The 
16' X 16' field has a pixel size of 0.24"; for each field four 
dithers were observed, each offset by 5" in order to estimate 
telescope shear estimation and to remove cosmetic defects. 
As for the Keck survey, we observed the fields as they passed 
through meridian to minimise flexure-induced distortions. 
Bias frames, sky flats and standard stars were acquired each 
night, and refocusing was carried out several times per night. 

The median seeing for the WHT fields is 0.8", which 
is adequate for our purposes (see Figure 8 of Bacon et al 
2001); Figure |^ records the seeing values for all fields. No 
data was used with seeing > 1.0". 

The dithers were exposed for 900 seconds each, amount- 
ing to a 1 hour exposure on each field. The median magni- 
tude of detected galaxies was i?=25.0, with an imcat signal- 
to-noise of 5.0 being reached for galaxies at i?=25.8, similar 
to our previous survey fields. We keep 15 resolved galax- 
ies per square arcminute in our final object catalogue (see 
section 3.4 for the selection criteria), which corresponds to 
a final magnitude median of 7? ~ 23.5 ± 0.2, including red- 
dening corrections. As with Keck, the distribution of magni- 
tudes for used galaxies extends substantially fainter than the 
median, with the galaxy count as a function of magnitude 
reaching 50% of maximum at J? = 24.4. Figure ^ records the 
reddening-corrected median magnitude vs number densities 
for all fields. According to e.g. Cohen et al (2000), the above 
median magnitude corresponds to a median source redshift 
of Zs ~ 0.8; using a similar estimate for the redshift error to 
that in section 3.1, we find an uncertainty in this quantity 
of 0.06. As this may be subject to sample variance, we will 
conservatively use a rounded uncertainty of 0.1. 



4.2 Data Reduction 

The reduction of the WHT fields followed the same stan- 



Figure 7. Example of astrometric distortions for both chips of 
the Wide Field Camera on WHT. 



subtraction and flat flelding were followed by the elimina- 
tion of fringing, which occurs on the flelds at 0.5% of sky 
background level. In order to defringe, all science exposures 
for each night were stacked without offsetting, using sigma- 
clipping to remove objects; this provided a fringe frame for 
each night. A multiple of the fringe frame which minimised 
rms background on each dither was subtracted, resulting in 
the fringes being reduced to an < 0.05% level, as in BRE. As- 
trometric matching and stacking of dithers then proceeded 
in the same manner as for the Keck data. 



4.3 Instrumental Distortions 

We re-checked the instrumental distortion for the WHT 



3.3 



dard procedure as outlined in BRE and §3.2 above. Bias 



flelds using the same method as for the Keck flelds (see 
and BRE). However, the larger number of objects in com- 
mon in the WHT exposures affords a determination of the 
instrumental shear on individual flelds rather than having 
to stack many fields as we did in 3.3, Figure |^ shows the 
WHT instrumental shear pattern for a typical field, includ- 
ing both chips of the new 16' x 16' mosaic field. Since the 
shear has a mean of 0.1% and is < 0.4% at the edges of the 
field, we again find that the correlation function of the tele- 
scope shear is negligible and need not be corrected for. The 
telescope shear estimates fluctuate by < 0.1% from fleld to 
fleld. Magnification and rotation components are also found 
to be < 0.4%, and are therefore negligible. 



4.4 Point Spread Function 

Correction for isotropic and anisotropic smear components 
proceeded in an entirely analogous fashion to that for the 
Keck data, using imcat 's ellipticities and polarisabilities to 
correct the galaxy shear estimates (see BRE for full details). 
The same criteria were used for removing noisy detections 
(rg > 1.0, > 15, e > 0.5). The rms ellipticity of our stars 
from fleld to field was (Te* — 0.05 (see BRE figure 7 for an 
example stellar ellipticity field). 

As with the Keck data, we checked for systematics in 
our shear estimators by taking all shear estimators from all 
fields and averaging them as a function of position (Figure 
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the initial selection cuts (see §3.5) are given equal weight 
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Figure 8. Average shear of all galaxies in all WHT fields as a 
function of position on chip. Overall, (71) = 0.00% it 0.13% and 
(72) =0.07% ±0.13%. 



|8|). For WHT we find that the mean shear for our entire en- 
semble is 71 = 0.0000±0.0013 and 72 = 0.0007±0.0013. This 
is entirely consistent with zero offset in the whole ensemble, 
as we would expect. Similarly, the figures demonstrate that 
there is no significant dependence of our shear values on 
position on the chip. 



5 COSMIC SHEAR STATISTICS: THE 
CORRELATION FUNCTION 

We are now in possession of good quality, wide area data 
for cosmic shear estimation from two telescopes. How are 
we to best extract the information contained from our shear 
estimator catalogues, in order to estimate cosmological pa- 
rameters and the amplitude of the mass power spectrum? 

In this paper we choose to use the shear correlation 
function (e.g. Kaiser 1998; Kamionkowski et al 1998), a use- 
ful probe of the shear power on various scales. It has the 
advantage of being simply related to the shear power spec- 
trum by a Fourier transform, and affords precise checks of 
the contribution to our results of PSF anisotropy systemat- 
ics. In this section, we shall describe the correlation functions 
and their use as estimators for the shear signal. 

The material which we have to work with is a set of cat- 
alogues containing the positions of galaxies, plus the shear 
estimator associated with each one; we therefore have a set 
of noisy samples of the cosmic shear field ^i{9). For this 
field, we can define shear correlation functions in the follow- 
ing way: 



C2{e) 



(7^(0)71(0)) 
(7^(0)72(0)), 



(1) 



where the average is to be taken over all galaxy pairs sepa- 
rated by an angle 9. In practice, of course, we will average 
over an annulus [6, 9 + AO]. Note that all objects that pass 



The shears used are those for a rotated coordinate system, 
where the a;-axis is defined by the line joining the two galax- 
ies. This can be achieved using the transformations 



7[ = 71 cos(2(^) + 72 sin(2</>) 
72 = -71 sin(20) + 72 cos(2<j!)). 



(2) 



where 7^ and 7J are the rotated shears that we require for 
our correlation functions. We can define a third correlation 
function 



C3 = {7r(o)75(e)) = {72(o)7r(e)), 



(3) 



which is expected to average to zero: under refiections C3 
changes sign, so we expect equal contributions from posi- 
tive and negative shear-shear configurations. Therefore C3 
provides a check on systematic effects introduced by our 
corrections, which need not be parity invariant. 

Finally, it is convenient to define a total correlation 
function 



C(0) = Ci{9) + C2{9). 



(4) 



The shear correlation functions can be readily calcu- 
lated for any cosmological model using 



C,{9) 



1 

47r 



dl IC^ 



(5) 



for i = 1,2, and where C7 is the shear power spectrum de- 
fined as in BRE (denoted as Ci in that paper). The correla- 
tion functions for a ACDM model with Qm = 0.3, F = 0.21, 
as = 1 and a median source redshift of Zs = 1.0 and 0.8 are 
plotted in figure ^ 

Since the fields are widely separated, we will not try 
to take the correlation functions for the whole area repre- 
sented. Instead, we measure the correlation functions in each 
separate field. However, we need to bear in mind that the 
measured correlation function in a given field is noisy due to 
the shot noise of the galaxy ellipticities. Thus we require an 
overall estimator of the correlation function of the cosmic 
shear for all fields. 

For this purpose, let us denote the correlation function 
measured in field / by C({6). An estimator for Ci{0) is sim- 
ply its average over all fields, i.e. 



Similarly, the uncertainty on our estimate of Ci{0) is simply 
the error in the mean given by 



(7) 



The values of Ci{9) at different are generally not inde- 
pendent. As a result, it is important for fair error estima- 
tion on cosmological parameters to have knowledge of the 
covariance of the correlation functions at different angular 
bins, cav[Ci{9),C'j{9')]. This can easily be calculated from 
the measured correlation functions using the relation 

coY[a{9),Cj{9')] ~ 
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Nf 

-T 



iie)-c,{e)] [cf{e')-c,ie')]. (8) 



The above expression has the advantage of containing 
information regarding the contribution to our errors of both 
shot noise and sample variance; that is to say, it contains a 
measure of the entire error budget for all scales, apart from 
systematic contributions. It also allows us to account for 
the covariance between different angular scales, as well as 
that between the different correlation functions. We can use 
these estimators for the correlation function, together with 
the covariance matrix, to find the best fit of cosmological 
parameters to our data. 

To assess the significance of the detection of the lensing 
signal, it is useful to consider the errors upon the correlation 
functions arising from statistical noise alone, i.e. neglecting 
sample variance. It is easy to show that the covariance ma- 
trix in this case is given by 



COV,tat[Q(e),C,(e') 



-^pairs 

(0) 



SijSggi, 



(9) 



where = (71) = {7I) is the intrinsic ellipticity variance 
of individual galaxies, and A^pairs(S) is the number of galaxy 
pairs used in the angular bin centered on 9. In using this 
equation, we will use the measured ellipticity dispersions of 
(jj ~ 0.31 for both WHT and Keck. Note that the covariance 
covstat vanishes for different correlation functions (i 7^ j) 
and for different angular scales {6 7^ 6'). 



5.1 Star-galaxy Correlation Functions 

We are able to apply the correlation function formalism to 
other quantities which we have measured besides the shear 
field. A particularly useful check of systematic effects is 
available by considering the extra contribution to the shear 
correlation function from PSF ellipticity contamination. If 
we have a small addition to the shear field due to uncor- 
rected contributions by the PSF ellipticity. 



7» = 7i + ae^ , 



(10) 



then it is clear that (7^6*) — a{e*e*) and {7i7i) = (jiyi) + 
(e*e*); from this it follows that the uncorrected elliptici- 
ties add a component to the measured correlation function 
described by 



cr = 



(7.er)' 



(11) 



where i = 1,2. This can be directly measured from our data 
to determine the error due to PSF systematics. 



5.2 Shear Variance 

Before applying these correlation functions to our data, we 
briefly examine their relationship to the cell-averaged shear 
variance; this is often used to quote and compare cosmic 
shear results, so it is important to be able to convert between 
the two estimators. 

Let us consider a cell with window function W{6) nor- 
malised as J d^e Wie) = 1. The average shear within this 
cell is 



7» 



d^0 wie)'y,{e). 



(12) 



It is easy to show that the shear variance within such a 
cell is related to the (total) shear correlation function by 



:i7i') 



(13) 



where we have used the small angle approximation. The 
shear variance is thus simply the average of the correlation 
function over the area of the aperture. It is also easy to show 
that the error variance in measuring cr^ is related to that for 
the correlation functions by 



2r 2i 



= d 



d^e' w{e)W{6')coY[C{e),c{e')]. (i4) 



For a circular cell of radius 6, we can calculate these 
integrals by noting that, for a top-hat window function, 
/ d^e'W{e') 261-2 ^Q,Q, ^^^^ ^j^ggg 

mulae to convert simply between correlation function and 
variance measurements. 



6 RESULTS 

6.1 Correlation Functions 

Now that we have equipped ourselves with the necessary 
tools for measurement, we proceed to examine the ampli- 
tude of the cosmic shear in our data. We first measured the 
correlation functions defined above, together with their co- 
variances, for all of our cosmic shear survey fields for Keck 
and WHT. Figure ^ compares the resulting shear correlation 
functions for both experiments, after removal of 3o" outliers 
in star-galaxy residual correlation. The inner error bars show 
the statistical errors derived from Equation (^). The outer 
error bars show the total error bar (shot noise -I- sample 
variance) derived from the diagonal elements of the covari- 
ance matrix in Equation (^) . Note that sample variance con- 
tributes significantly to our uncertainties and should there- 
fore be included when constraining cosmological parameters. 
The off-diagonal elements (see Figures ^ and |l^) allow us to 
quantify the correlation between the different angular bins 
and between the different correlation functions. 

The expected correlation function in a ACDM model 
with Q,m ~ 0.3 is also shown in Figure ^ The shape pa- 
rameter was set to F = 0.21, close to the values indicated 
by recent galaxy surveys (Percival et al. 2001; Szalay et al. 
2001) while maintaining the F ~ ^rnh relation with h = 0.7. 
The theoretical curves were normalised by taking as = 1, 
consistent with 'old' cluster normalisations (Eke et al. 1998; 
Viana & Liddle 1999; Pierpaoh et al. 2001). A comparison 
with the newer cluster normalisation of erg = 0.7 (Borgani 
et al 2001; Seljak 2001; Reiprich & Bohringer 2001; Viana 
et al. 2001) will be presented below. The median redshift for 
the galaxies was assumed to be Zs = 0.8 and 1.0, as relevant 
for WHT and Keck respectively (see ^ and §^). For the 
models, the redshift distribution of the galaxies was taken 
to be as in BRE. 

Importantly, despite different strategies and instru- 
ments, the WHT and Keck measurements of Ci and C2 are 
in good agreement, especially after taking account of the 
difference in z^. Moreover, both agree with the displayed 
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Figure 9. Shear correlation functions compared for the WHT 
and Keck. The inner error bars correspond to statistical errors 
only, while the outer error bars correspond to the total errors 
(statistical + sample variance). For comparison, the correlation 
functions expected for a ACDM model with H™, = 0.3, erg = 1 
and r = 0.21. A median source redshift of 2^=0.8 and 1.0 are 
shown, as appropriate for WHT and Keck, respectively. 



ACDM models (see below for detailed discussion) . Our mea- 
surement of C3 is shown on the bottom panel and is consis- 
tent with zero as expected from parity s ymm etry. 

We can use the formalism of section 5.1 to evaluate the 



systematic contribution of residual stellar anisotropy to our 
shear correlation function. For this purpose, we measure the 
star-galaxy correlation function given by equation (^l|) for 
all of our WHT fields, and show the results on Figure |lo[ 
Note that, bearing in mind that the scale on this figure is 
the same as for our shear correlation function, the system- 
atic contribution due to our anisotropy correction is entirely 
negligible; it is lower than the shear signal by a factor > 10 
everywhere. 



6.2 Sources of Error and Covariance between 
Angular Scales 

It is instructive to compare the different sources of error 
for each telescope, given the complementary survey strate- 
gies. Figure |ll| shows the error variance of the shear correla- 
tion correlation functions Ci -I- C2 of Figure ^ as a function 
of angular scale d. The total error variance (from Eq. Q) 
was decomposed into a statistical error (from Eq. |^) and 
a sample variance contribution (computed by subtracting 
the latter from the former). The cosmic variance contribu- 
tion is clearly significant for both Keck and WHT and must 
be taken into account in the determination of cosmological 
parameters from cosmic shear data. 

The advantage of using a lOm-class telescope in reduc- 
ing statistical errors is apparent in this Figure. Indeed, the 



Figure 10. Star-galaxy correlation contribution to the error 
budget. Demonstrates the negligbly small residual alignment of 
galaxy shear measurements to the orientation of the original PSF 
(note the same scale on the j/-axis as figure ^ above). 
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Figure 11. Errors in the shear correlation correlation function 
C\ -f C2 as a function of angular scale 9 for both WHT (top) and 
Keck (bottom) . The total error variance (solid line) is decomposed 
into a statistical (dotted line) and a sample variance (dashed line) 
contribution in each case. Note that the area of the WHT survey 
is a factor of about 1.7 larger than that for Keck. 
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Figure 12. Covariance matrix of the shear C\ correlation func- 
tion between the different angular bins of the WHT data plotted 
in figure |9| 

statistical errors for Keck for ;^ 3' are significantly smaller 
than that for WHT, after scaling them down by a factor of 
1.7 corresponding to a common survey area of 1 deg^ in both 
surveys. The much shorter exposure times allowed by Keck 
yield an improvement in the seeing, tracking errors and over- 
all image quality. Aided by the finer pixellisation of the ESI 
CCD, this results in an increased galaxy number density and 
correspondingly lower statistical errors when normalised by 
survey area (as predicted by Bacon et al 2001). Of course, 
the statistical errors for WHT are smaller on large angu- 
lar scales [6 ^ 3') thanks to the larger field of view of this 
telescope. 

The covariance of the shear correlation functions on dif- 
ferent angular scales is also non-negligible. Figures |l|and|l| 
show the covariance matrix for the Ci correlation function 
for WHT and Keck, respectively. In each case, the angu- 
lar bins correspond to the 9 values shown in Figure The 
WHT correlation function clearly has important covariance 
between bins on small angular scales (6 ^ 3'). The Keck 
correlation functions also has important covariance on large 
scales {9 ^ 3'). This is due to the elongated field geome- 
try of the Keck camera. In both cases, it is important to 
include the full covariance matrix for the determination of 
cosmological parameters. 

6.3 Shear Variance 

We now wish to compare our results to those obtained re- 
cently by other groups. Since many of these results are 
quoted in terms of the angular dependence of the shear 
variance, we use equation ( |l3| ) to convert our correlation 
function measurements into shear variance measures as a 
function of scale. 




0.0 0.5 1.0 1.5 2.0 



COV[C.(0).C.(0O] 




e bin 

Figure 13. Covariance matrix of the shear C'l correlation func- 
tion between the different angular bins of the Keck data plotted 
in figure ^ 

The results for the WHT and Keck surveys, are plot- 
ted in Figure |l^. Taking the most significant point in each 
case and using statistical errors only, we find that the cosmic 
shear signal is detected at the 3.7a, 3.5a and 5.1a level with 
Keck, WHT, and both combined, respectively. Also plotted 
are the predictions for the ACDM model with Qrn ~ 0.3 
and r — 0.21 as before. They are shown for a range of val- 
ues of the galaxy median redshift Zs , corresponding roughly 
to the uncertainty and survey-to-survey variations in this 
parameter. Both old [ag = 1) and new cluster normali- 
sations {as ~ 0.7) are shown (see discussion below). We 
also plotted the measurements from other groups, namely 
from van Waerbeke et al 2000 (CFHT vW-|-), Kaiser et al. 

2000 (CFHT K+), Wittman et al 2000 (CTIO), MaoU et al 

2001 (VLT), and van Waerbeke et al 2001 (CFHT vW4-+). 
Note that the errors for the VLT, CFHT vW-l- and CFHT 
vW-l— I- do not include cosmic variance. All of the above 
groups are obtaining results which are broadly consistent 
with each other and with ACDM with the older normalisa- 
tion of as = 1. 



6.4 Cosmological Constraints 

Now that we have measurements of cosmic shear on sev- 
eral scales, we can determine the implications for cosmo- 
logical parameters. In order to do this, we use a Maxi- 
mum Likelihood approach. We first construct data vectors 
d = {C\(9n),C-i.(9n)\ which are simply a rearrangement of 
our observed correlation functions. We consider a ACDM 
model with two parameters, ai = f2m and 02 = as. The 
shape parameter F for the matter power spectrum is set to 
0.21 as indicated by recent measurements of galaxy cluster- 
ing (Percival et al. 2001; Szalay et al. 2001). Using the for- 



© 2002 RAS, MNRAS 000, 



Joint Cosmic Shear Measurements with the Keck and William Herschel Telescopes 11 



0.0008 



0.0006 



0.0004 



0.0002 



0.0000 



z,= 1.0,0.9,0.8 
ACDM, a,= l.Q 
ACDM, C7.= 0.7 



CFHT vWh 
X CFHT vWh 
C CTIO 

CFHT K+ 
~ VLT 

awht 

□ Keck 




1.5 



. 1.0 



0.5 - 





1 1 1 


1 1 1 1 1 1 1 
- Keck 








- Clusters 


fold) 


■ ^1 v 


\\\ 


■■■ Clusters 


(new; 










■ , , , 


. . 1 


. . 1 , . . 1 


, , 1 , , ,■ 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 14. Shear variance as a function of aperture radius. Our 
measurements with Keck and WHT are shown along with those 
of other groups (see text). Again, the inner and outer error bars 
correspond to statistical and total uncertainties, respectively. The 
expected shear variances for a ACDM model with Qrn = 0.3, and 
r = 0.21 are also shown for a range of values for the median source 
redshift Zs corresponding approximately to the uncertainty in this 
parameter. The models are shown both for erg ~ 1 (solid lines) 
and ag ~ 0.7 (dashed) cluster normalisation. 



Figure 15. Constraints on the joint distribution of Qrn and erg 
from our Keck measurement. The solid contours correspond to 
68% and 95% confidence levels. Note that these contours include 
statistical errors and non-gaussian sample variance, but do not 
include the uncertainty in the galaxy redshifts and in the shear 
measurement method. Also shown are the Itr contours from the 
old (Pierpaoli et al. 2001) and new (Seljak 2001) cluster abun- 
dance normalisation. 



malism described in BRE and in section g, we can compute 
the correlation function for any values of these parameters 
and arrange them as a theory vector t(a) set out in exactly 
the same form as our data vector. The median redshift Zs is 
set to 0.8 and 1.0 for WHT and Keck respectively. 

Because the correlation functions were derived from an 
average over a large number of fields, the central limit theo- 
rem ensures that our errors (and covariances) are gaussian. 
In this case the log-likelihood is simply 

x' = [d-t(a)]^V-^[d-t(a)], (15) 

where V is the covariance matrix computed using Equa- 
tion (^). We minimise this quantity as a function of a to 
find the best estimate of the normalisation of the power spec- 
trum. To compute confidence contours, we numerically inte- 
grate the probability distribution function P(a) cx 

The constraints for our Keck and WHT data, taken sep- 
arately, are shown on figures ^ and |l6| Note that they do 
not include the uncertainty in Zs and in the shear measure- 
ment method. These two sources of errors will be included 
in our final estimate below. The constraints from each of the 
data set are consistent with each other. As is apparent on 
the figures, the constraints reveal the well known degeneracy 
between as and ilm when 2-point statistics alone are used. 
The wider angular range of WHT allow us to start break- 
ing this degeneracy, by rejecting large values of Qm (see van 
Waerbeke et al. 2001). Note that the width of the contours 
for Keck and WHT are comparable, even though the WHT 
area is larger by about 70%. This results from the larger 
number density of galaxies afforded by Keck thanks to the 
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Figure 16. Constraints on the joint distribution of Qrn and rrg 
as in the previous figure, but for our WHT measurement. 



better seeing statistics and smaller pixel size (see discussion 
in Bacon et al. 2001). 

The constraints on ag and obtained by combining 
the Keck and WHT data are shown in Figure ^ As ex- 
pected, the combined contours are consistent with the two 
measurements taken separately. A good fit to our 68% con- 
fidence level is given by 
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Figure 17. Constraints on the joint distribution of f2m and erg 
for the combination of our Keck and WHT measurements. 




= 0.97_o'o9 (shot noise + sample variance)(16) 



with 0.14 < Q,m < 0.65. To this statistical error (which 
includes non-gaussian sample variance), we also must add 
that arising from the uncertainty the in the redshift dis- 
tribution and in the method calibration. A good estimate 
for these uncertainties can be obtained by noting that the 
shear variance on the central scale for our experiments, i.e. 
circular radius of 3 arcmins, scales as ct^ oc P~'^ag^zl'^, 
where P^^ is the ellipticity-to-shear conversion factor in the 
KSB method (see BRE). This factor was shown to have a 
la uncertainty of 5% (BRE) . The uncertainty in Zs is about 
12% for our two data sets (see §^ and §^). Using the for- 
mer expression to propagate the errors yields a final con- 
straint on the amplitude of the matter power spectrum of 
cr8(n™/0.3)°-^* = 0.97to.M ± 0.07 ± 0.04, where the error 
bars denote 68% confidence levels, arising from statistical 
(incl. sample variance), redshift uncertainty and uncer- 
tainty, respectively. Adding these errors in quadrature, we 
obtain our final estimate for the amplitude of the power 
spectrum, namely 

/O \0-6S 

(78 f-^ J =0.97 ±0.13 (total error), (17) 

where the error is the total error for a 68% confidence level. 

Our results are consistent with the results of van 
Waerbeke et al. (2002) who found erg = (0.94 ± 
0.07)(n^/0.3)°-^''^o « (68%CL), after marginalising 

over r in the range 0.1 to 0.4. They are also consistent with 
those of Hoekstra et al. (2002) who found ag = 0.8ll[5:i9 
(95%CL) for = 0.3 and F = 0.21. 

It is interesting to compare our measurement of erg 
to that derived from cluster abundance. Our results are 
in agreement with the older normalisation of the power 
spectrum from cluster abundance. Until recently, the clus- 
ter abundance was indeed quoted to be erg = (1.00 ± 
0.07)(0.3/!:!„0°'^''(r/0.21)''-28-0-3i"™ (la; Pierpaoh et al. 



2001; see also earlier and consistent estimates by Eke et al. 
1998; Viana and Liddle 1999). Recently the cluster normal- 
isation has been revised to a lower value, mainly because 
of the use of the observed mass-temperature relation for 
clusters, rather than the simulated one. For instance, Sel- 
jak (2001) finds as = (0.75 ± 0.06)(0.3/n™)° **(F/0.21)° °* 
(Icr; see also Borgani et al 2001; Rieprich & Bohringer 2001; 
Viana et al. 2001). A similar normalisation was found by 
combining the 2dF galaxy survey with CMB measurements 
by Lahav et al. (2001). The two cluster normalisations are 
also plotted on Figures |l5[ |l| and 0. Clearly our results 
are most consistent with a trg = 1 normalisation, but do not 
rule out the new normalisation which is consistent at the 
1.7a level for ^Im = 0.3. This can also be seen by exam- 
ining Figure |l^ which shows the impact of these different 
normalisations for the shear variance. 



7 CONCLUSIONS 

We have presented a thorough analysis of weak lensing sig- 
nals obtained using two telescopes. The independent in- 
struments introduce different systematic effects upon galaxy 
shapes, the cross-checking and control of which is the most 
important challenge facing current measurements of weak 
shear. Our many-line-of-sight survey strategy for cosmic 
shear, using 173 Keck telescope fields, also complements an 
alternative 20 large-field survey strategy with the William 
Herschel Telescope. 

We have measured cosmic shear at a signal-to-noise of 
5.1a with both Keck and WHT, and have measured the am- 
plitude of the power spectrum, erg (§3") = 0.97 ± 0.13 
with 0.14 < rim < 0.65, including all contributions to the 
68% CL uncertainty: statistical noise, sample variance, co- 
variance between angular bins, systematic effects, and red- 
shift uncertainty. A measurement of this quantity from cos- 
mic shear is cosmologically valuable, as it represents a direct 
measure of the amplitude of mass fluctuations. 

These measurements have been obtained after a care- 
ful study and removal of systematic effects which can mimic 
a shear signal at the 1% level. We have demonstrated for 
both telescopes that no offset in the ensemble of shear es- 
timators is found as a function of position on images, and 
that the contribution of star-galaxy correlations is negligible 
with our catalogue selection. Our methods have been tested 
thoroughly in the context of detailed simulations of realistic 
images (Bacon et al. 2001) and have been shown to operate 
successfully in recovering shear at the necessary level. 

Our results for Keck and WHT are consistent with each 
other, strengthening confidence in control of systematics. 
The joint results are also consistent with other recent cos- 
mic shear measurements (Hoekstra et al. 2002; van Waer- 
beke et al. 2002) . They also agree with the old cluster abun- 
dance normalisation (Eke et al. 1998; Viana and Liddle 1999; 
Pierpaoli et al. 2001). Our results prefer this normalisation, 
but can not rule out lower cluster-abundance normalisation 
which has recently been derived (Borgani et al 2001; Seljak 
2001; Rieprich & Bohringer 200; Viana et al. 2001; see also 
the similar normalisation derived from 2dF-|-CMB by Lahav 
et al. 2001). This discrepancy, if confirmed, could arise from 
unknown systematics in either the cluster or cosmic shear 
methods. Note for instance that uncorrected systematics in 
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cosmic shear measurements will tend to add to the lens- 
ing signal and thus lead to an overestimation of ag. For the 
cluster method, further studies would be needed to under- 
stand the difference between the observed mass-temperature 
relation and that found in numerical simulations. It is im- 
portant to understand the origin of the discrepancy between 
cosmic shear and cluster abundance methods. If this is not 
explained by such systematics, it could point towards a fail- 
ure of the standard ACDM paradigm, and therefore have 
important consequences for cosmology. 
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